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Abstract 

The  Naval  Surface  Warfare  Center,  Dahlgren 
Division  Aeroprediction  Code  has  been  extended  to 
the  roll  position  of  45  deg  (or  fins  in  cross  or  “x” 
orientation).  New  technology  developed  and 
discussed  in  this  paper  includes:  wing-body  and 
body-wing  interference  factors  for  this  roll  orientation 
and  to  angles  of  attack  of  90  deg;  approximate 
methods  to  estimate  wing-alone  center  of  pressure 
shift  at  high  angles  of  attack  and  45  deg  roll  position; 
and  “fin  choking”  at  high  Mach  number.  This  new 
technology  allows  aerodynamics  to  be  estimated  with 
average  accuracy  levels  of  ±10  percent  on  normal  and 
axial  force  coefficients  and  ±4  percent  of  body  length 
on  center  of  pressure.  Exceptions  to  this  are  at 
subsonic  Mach  number  and  high  angle  of  attack 
where  wind  tunnel  sting  interference  effects  are 
present,  and  at  high  Mach  number  and  angle  of  attack 
where  internal  shock  interactions  from  a  forward  fin 
onto  an  aft-mounted  fin  are  present.  Results  are 
presented  for  several  wing-body  and  wing-body-tail 
configurations  at  various  Mach  numbers  and  angles 
of  attack  to  support  this  average  accuracy  level 
conclusion. 


Nomenclature 

AR  Aspect  Ratio  =  b2/Aw 

a  Speed  of  Sound  (ft/sec) 

Aref  Reference  area  (maximum  cross- 

sectional  area  of  body,  if  a  body  is 
present,  or  planform  area  of  wing, 
if  wing  aloneXft2) 

Aw  Planform  area  of  wing  in  crossflow 

plane  (ft2) 

b  Wing  span  (not  including  body)(ft) 

C .  Axial  force  coefficient 

A 
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CM 

Pitching  moment  coefficient  (based 
on  reference  area  and  body 
diameter,  if  body  present,  or  mean 
aerodynamic  chord,  if  wing  alone) 

cN 

Normal  force  coefficient 

CN„ 

Normal  force  coefficient  of  body 
alone 

^nbcv> 

Negative  afterbody  normal-force 
coefficient  due  to  canard  or  wing- 
shed  vortices 

Cy  »  CN 
nB(W>  ”b<t> 

Normal-force  coefficient  on  body 
in  presence  of  wing  or  tail 

ACNB(W, 

Additional  normal-force  coefficient 
on  body  due  to  presence  of  wing 

C»L 

Linear  component  of  normal-force 
coefficient 

CNnl 

Nonlinear  component  of  normal- 
force  coefficient 

CN 

WTCV) 

Negative  normal-force  coefficient 
component  on  tail  due  to  wing  or 
canard-shed  vortex 

CNW 

Normal  force  coefficient  of  wing 
alone 

CN 

^W(B) 

Normal-force  coefficient  of  wing  in 
presence  of  body 

CN« 

Normal-force  coefficient  derivative 

Cr 

Root  chord  (ft) 

c, 

Tip  chord  (ft) 

c„  c2 

Dimensionless  empirical  factors 
used  in  nonlinear  model  of  kW(B)  to 
approximate  effects  due  to  high 
AOA  or  control  deflection 

cal 

Caliber(s)  (one  body  diameter) 

DUO  QUALHs 


UNCLASSIFIED 


9961213  005 


UNCLASSIFIED 


d,ef 

Reference  body  diameter  (ft) 

Kb(w)?  K-b(t) 

Ratio  of  additional  body  normal- 
force  coefficient  in  presence  of 
wing  or  tail  to  that  of  wing  or  tail 
alone  at  8  =  0  deg 

^B(W)5  kB(T) 

Ratio  of  additional  body  normal- 
force  coefficient  due  to  presence  of 
wing  or  tail  at  a  control  deflection 
to  that  of  wing  or  tail  alone  at 
oc  =  o  deg 

[KB(W)]mIN 

Minimum  value  of  KB(W)  as  percent 
of  slender-body  theory  value 

Kb, 

Lowest  value  of  [K^]^  for 
=  45  deg 

Kw(B),  ^t<b) 

Ratio  of  normal-force  coefficient  of 
wing  or  tail  in  presence  of  body  to 
that  of  wing  or  tail  alone  at  8  =  0 
deg 

k\V(B)>  LT(B) 

Ratio  of  wing  or  tail  normal-force 
coefficient  in  presence  of  body  due 
to  a  control  deflection  to  that  of 
wing  or  tail  alone  at  «  =  0  deg 

Length  of  body  (cal  or  ft) 

Afterbody  length  (cal  or  ft) 

M 

Mach  number  =  V/a 

r 

Radius  of  body  (ft) 

s 

Wing  or  tail  semispan  plus  body 
radius  in  wing-body  lift 
methodology 

V 

Velocity  (ft/sec) 

XCP 

Center  of  pressure  (in  feet  or 
calibers  from  some  reference  point 
that  can  be  specified)  in  x  direction 

oc 

Angle  of  attack  (deg) 

“w=  “T 

Local  angle  of  attack  of  wing  or  tail 
(«  +  8  in  degrees) 

6 

Control  deflection  (deg),  positive 
leading  edge  up 

X 

Taper  ratio  of  a  lifting  surface  = 
c/cr 

<t>  Circumferential  position  around 

body  where  5>  =  0  is  leeward  plane 
with  fins  in  plus  fm  arrangement 

(deg) 

oo  Free-stream  conditions 

Introduction  and  Background 

Many  of  the  world’s  missiles  fly  in  either  the  roll 
stabilized  position  of  <X>  =  0  deg  (or  plus  fin 
orientation)  or  O  =  45  deg  (or  cross  fin  orientation). 
Figure  1  illustrates  these  fm  orientations  for  a 
cruciform  missile  looking  from  the  missile  nose 
toward  the  rear.  As  illustrated  in  the  figure,  the 
0  =  0  deg  plane  generally  gives  slightly  more  normal 
force  and  a  slightly  more  stable  configuration  in  pitch 
at  a  given  angle  of  attack  (AOA)  than  does  the  missile 
rolled  to  O  =  45  deg  (the  physics  of  why  this  occurs 
will  be  discussed  in  the  analysis  section).  On  the 
other  hand,  a  missile  in  the  O  =  45  deg  plane  is  in  a 
roll-stable  position,  which  means  less  energy  is 
required  to  maintain  a  constant  roll  orientation.  Also, 
all  four  fins  can  be  deflected  simultaneously,  giving 
30  to  50  percent  more  normal  force  from  control 
deflection  than  only  two  fins  deflected  in  the  O  =  0 
deg  roll  position. 

The  latest  version  (1995)1  of  the  Naval  Surface 
Warfare  Center  Dahlgren  Division  (NSWCDD) 
Aeroprediction  Code,  AP95,  calculates  aerodynamics 
only  in  the  O  =  0  deg  plane.  As  such,  aerodynamics 
at  O  =  45  deg  roll  must  be  obtained  from  another 
aerodynamics  code2,3  or  estimated  from  the  AP95 
0  =  0  deg  results.  It  is  the  intent  of  this  paper  to 
discuss  the  physical  phenomena  that  occur  at  the 
0  =  0  and  45  deg  roll  orientations  and  to  summarize  a 
semiempirical  mathematical  model  that  allows  the 
AP95  to  be  extended  to  the  O  =  45  deg  plane. 

Details  of  this  work  can  be  found  in  Ref.  4. 

The  AP95  is  an  approximate  analytical 
aeroprediction  code  primarily  designed  to  provide 
preliminary  estimates  of  aerodynamics  for  use  in 
particle  ballistic  models,  trim  aerodynamic  models,  or 
structures  and  heat  transfer  models.  Static 
aerodynamics  are  generally  estimated  with  average 
accuracy  levels  of  ±10  percent  and  center  of  pressure 
(COP)  within  ±4  percent  of  the  body  length.  To 
obtain  aerodynamics  estimation  accuracy  levels 
generally  desired  for  full  six-degree-of-freedom 
simulations  requires  either  a  more  accurate  numerical 
code5  or  wind  tunnel  data  or  both.  As  such,  the  AP95 
does  not  attempt  to  compute  out-of-the-pitch-plane 
aerodynamics  or  coupling  effects  between  the  pitch 
and  yaw  planes.  On  the  other  hand,  if  one  is 
interested  in  the  =  45  deg  plane,  and  control 
deflections  are  symmetric  with  respect  to  the  pitch 
plane,  then  the  AP95  code  can  be  modified  to  allow 
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0=0°  0=45° 


|  VooSin  a 


|  VMsin  a 


TYPICAL  FORCE  AND  MOMENT  COMPARISONS 


0=0! 


0=45° 


•  (Cn)<d=oo  generally  HIGHER  THAN 
(Cn)<i>=45°  AT  HIGHER  a 

*  GENERALLY  MORE  STABLE 
IN  PITCH  AT  0=0°  vs  0=45° 


•  NATURALLY  STABLE  POSITION 
IN  ROLL  (LESS  ENERGY  TO 
MAINTAIN  CONSTANT  ROLL) 

•  MORE  CONTROL  AUTHORITY 
AT  0=45°  DUE  TO  FOUR  FINS 
DEFLECTED  vs  TWO  AT  0=0° 


FIGURE  1.  CRUCIFORM  WING-BODY-TAIL  MISSILE  CONFIGURATION 
FLYING  AT  ROLL  OF  O  =  0  DEG  AND  45  DEG 
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aerodynamics  for  the  <£  =  45  deg  plane.  These 
modifications  will  be  made  in  an  analogous  process 
to  the  nonlinear  semiempirical  methodology  of 
Ref  L 

The  overall  approach  to  modify  the  AP95  code  for 
nonlinear  aerodynamics  in  the  O  =  45  deg  plane  will 
thus  be  very  similar  to  that  for  the  O  =  0  deg  plane.1,6 
Linearized  theories  or  slender  body  theory  (SBT)  will 
be  used  for  low  AOA  estimates  and  the  data  bases  of 
Refs.  7-10  will  be  used  to  develop  empirical  or 
semiempirical  corrections  to  account  for  the 
nonlinearities  that  occur  in  normal  force  and  COP 
with  increasing  AOA. 


Analysis 

As  indicated  in  the  introduction,  the  goal  is  to 
develop  a  nonlinear  semiempirical  model  for 
cruciform  missiles  for  the  3>  =  45  deg  plane.  It  is 
envisioned  that  this  model  will  be  analogous  to  the 
0  =  0  deg  plane  methods  in  AP95,  except  the  normal 
force  and  pitching  moments  due  to  the  wing  alone 
and  interference  aerodynamics  will  have  to  be  derived 
for  the  O  =  45  deg  roll  orientation. 

Referring  to  the  total  normal  force  coefficient 
equation  for  a  wing-body-tail  configuration  as  given 
by  Ref  10,  we  have 

CN  =  CNb  +  [(KW(B)  +  KB(W))  06 
+  fcw(B)  +  kB(W))  5W  ]  (CNoc)w 

+  [(^T(B)  +  KB(T))  **  +  ^T(B)  +  kB(T))  ^t]  (CN<x)t 


The  CN  term  of  Eq.  (1),  which  is  the  downwash 

normal  force  on  the  body  due  to  the  wing  shed 
vortices,  is  neglected.  This  is  because  it  is  inherently 
included  in  the  wind  tunnel  data  bases,  and  it  is 
believed  the  errors  in  trying  to  analytically  estimate 
the  term,  subtract  it  out  on  one  configuration,  and 
then  add  it  back  in  later  on  a  different  configuration, 
are  as  large  or  larger  than  the  errors  from 
incorporating  it  into  the  KB(W)  term. 

Equation  (1)  can  also  be  rewritten  as 


CN  =  CN  +  CN  +  CN 

N  NB  WW(B)  ^BfW) 


+  cN  +  cN  +  cN 

WTCB>  ^BfT)  WT(V) 


(2A) 


where  it  is  understood  that  C*.  encompasses  the 

WBCT) 

r  term.  For  ease  of  implementation  into  an 

NB(V) 

existing  code  designed  primarily  for  linear 
aerodynamics,  most  of  the  terms  in  Eq.  (2A)  are 
separated  into  a  linear  and  nonlinear  contribution  due 
to  «  or  8.  For  example,  the  wing-body  term  is 
computed  in  the  AP95  code  as  follows: 


CNw(B)  [(CnJL  +  (CnJNL  w 


The  linear  or  small  AOA  terms  of  Eq.  (2B)  are 
estimated  by  linear  theory  (LT)  or  SBT.  This  gives 
the  Aeroprediction  Code  a  good  fundamental  basis 
for  its  aerodynamic  estimates.  The  nonlinear 
corrections  due  to  higher  AOA  or  control  deflection 
are  each  estimated  directly  from  component  wind 
tunnel  data  bases.7'10  Each  of  the  other  terms  in 
Eq.  (2A)  is  treated  in  a  similar  fashion  to  Eq.  (2B)  in 
the  actual  implementation  into  the  Aeroprediction 
Code. 

In  the  context  of  Eq.  (2),  we  therefore  seek  the 
nonlinear  definition  of  each  of  the  terms  in  Eq.  (2)  for 
=  45  deg  roll.  It  is  expected  that  the  body-alone 
term  [first  term  of  Eq.  (2  A)]  will  be  independent  of 
<J>.  In  reality,  this  is  not  necessarily  the  case  for 
M  <  2  and  high  AOA  because  of  the  asymmetric 
shedding  of  vortices.  The  mechanism  of  this 
shedding  is  not  clear,  but  it  is  suspected  that  slight 
imperfections  in  the  flow  or  body  shape,  from 
uniform  or  axisymmetric,  respectively,  could 
contribute  to  this  phenomenon.  At  present,  the 
Aeroprediction  Code  does  not  account  for  out-of- 
plane  aerodynamics,  and  therefore  the  side  force 
created  by  the  asymmetric  shedding  of  body  vortices 
is  not  predicted.  Also,  in  the  Ref.  9  data,  normal 
force  varied  by  about  10  percent  as  a  function  of  roll 
in  the  region  of  asymmetric  vortex  shedding.  Instead 
of  including  this  variation,  it  was  averaged  out. 

Each  of  the  remaining  terms  in  Eq.  (2)  will  be 
predicted  in  an  analogous  fashion  to  the  AP95 
developed  for  O  =  0  deg,  except  here  the  quantities 
will  be  for  <D  =  45  deg.  As  already  mentioned,  that 
approach  was  based  on  LT  or  SBT  for  small  values  of 
«  and  empirical  data  bases  (Refs.  7-10)  to  develop 
nonlinear  corrections  for  large  <*.  As  such,  it  is 
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instructive  to  examine  the  fundamental  impact  of  roll 
orientation  on  LT  and  SBT  before  proceeding  to  the 
nonlinear  corrections,  which  are  empirical  in  nature. 

SBT  and  LT  Results  for  Roll-Dependent 
Aerodynamics 

References  1 1  and  12  were  primary  materials  used 
for  examining  roll  dependence  implications  from 
slender  body  and  linearized  theories.  A  somewhat 
detailed  summary  of  these  results  is  given  for 
information  purposes  in  Ref.  4.  The  summary  of  the 
key  findings  in  Ref.  4,  repeated  here  for  convenience, 
are: 

a)  For  cruciform  wings  alone  or  a  wing-body 
combination,  the  total  normal  force  is 
independent  of  roll. 

b)  For  a  planar  wing-body  combination  at  roll, 
the  loading  on  the  windward  plane  panel  is 
greater  by  an  amount  equal  to  that  on  the 
leeward  plane  panel.  This  means  that  if  one 
were  trying  to  design  a  code  for  lateral 
aerodynamics,  roll  dependence  of  each  fin 
planform  must  be  considered.  On  the  other 
hand,  if  longitudinal  aerodynamics  are  of 
primary  interest,  the  total  normal  force  on 
the  entire  wing  planform  can  be  considered. 

c)  For  a  cruciform  wing-body-tail  configuration 
at  roll,  eight  vortices  are  shed  in  the  wing- 
body  region,  which  adversely  affects  the  tail 
lift.  This  is  as  opposed  to  four  vortices  at 

0  =  0  deg. 

d)  The  planar  theory  developed  for  wing-tail 
interference  can  be  used  to  approximate  the 
loss  of  lift  on  the  tails  at  O  =  45  deg. 

e)  The  aerodynamics  of  a  cruciform  wing- 
body-tail  combination  with  zero  control 
deflections  are  independent  of  roll  position. 

These  findings  for  roll  dependence  from  LT  or 
SBT  are  quite  useful  in  helping  plan  how  to  develop  a 
nonlinear  aeroprediction  code  for  O  =  45  deg.  While 
the  conclusions  of  LT  roll  dependence  may  not 
translate  to  the  nonlinear  case,  we  will  still  use  the 
findings  to  help  guide  the  nonlinear  code 
development.  In  particular,  the  item  (a)  conclusion 
implies  use  of  the  0  =  0  deg,  wing-alone  data  for 
O  =  45  deg.  This  is  quite  important  because  the 
available  wing-alone  data  bases  are  all  at  O  =  0  deg. 
This  means  that  any  nonlinear  wing-alone  roll 
dependence  will  be  included  in  the  interference 
factors  rather  than  the  wing-alone  solution,  which  is 
independent  of  O. 


The  second  major  result  of  the  key  SBT/LT  roll 
dependence  findings  is  that  for  cruciform  missiles,  we 
can  use  the  same  interference  approaches  as  in  the 
AP95,  except  the  constants  need  to  be  changed 
because  of  a  different  roll  angle.  The  combination  of 
these  two  conclusions  are  quite  important  because 
they  basically  allow  the  direct  usage  of  the  AP95 
code  with  different  constants  for  the  nonlinear 
interference  terms  at  O  =  45  deg  versus  0  =  0  deg. 

The  third  significant  conclusion  is  that  for  small 
AOA,  wing-body-tail  aerodynamics  are  independent 
of  roll  position.  This  allows  the  usage  of  wing-tail 
interference  methodology  designed  for  planar 
computations  for  different  roll  orientations,  so  long  as 
the  proper  number  of  vortices  are  considered.  Again, 
different  nonlinear  corrections  are  expected  for  the 
O  =  45  deg  versus  the  O  =  0  deg  roll  position. 

Nonlinear  Aerodynamics  Methods 

This  section  will  describe  the  methods  used  for 
computing  the  nonlinear  corrections  for  each  of  the 
terms  in  Eq.  (2).  These  corrections,  with  the 
exception  of  the  body  alone,  are  all  empirical  in 
nature. 

The  wing-  and  body-alone  methods  for  normal 
force  coefficient  prediction  are  quite  similar  to  those 
in  Ref.  1.  In  Ref.  1,  the  wing-alone  normal  force  was 
estimated  by  a  fourth  order  equation  in  AOA,  with  the 
constants  chosen  by  the  data  bases  of  Refs.  2,  7,  and 
8.  The  body-alone  normal  force  was  predicted  by 
linearized  theory  combined  with  a  modified  version 
of  the  Allen  Perkins  viscous  crossflow  theory13  at 
higher  AOA. 

The  COP  of  the  wing  at  0  =  0  deg  roll  was 
estimated  by  LT  at  low  AOA  and  assuming  the  COP 
goes  to  the  centroid  of  presented  area  at  AOA  60  deg. 
The  body-alone  COP  at  =  0  deg  was  estimated  by  a 
weighted  average  of  linearized  theory  and  viscous 
crossflow  theory  with  COP  shifts  based  on  data  at 
transonic  Mach  numbers.  The  3>  =  45  deg  body- 
alone  COP  is  assumed  to  be  the  same  as  that  at 
$  =  0  deg. 

A  shift  in  the  wing-body  COP  at  $  =  45  deg  has 
been  derived.4  This  shift  is  driven  by  the  asymmetric 
loading  that  occurs  on  the  windward  to  leeward  plane 
fins  as  AOA  is  increased. 

To  visualize  this  effect,  imagine  a  missile  rolled  to 
<I>  =  45  deg  and  increasing  in  AOA.  As  AOA 
increases,  two  things  occur.  First,  the  windward 
plane  fins  carry  more  and  more  of  the  load  compared 
to  the  leeward  plane  fins.  Second,  the  local  Mach 
number  in  the  windward  plane  is  different  and 
typically  lower  than  the  leeward  plane.  This  has  the 
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effect  of  shifting  the  wing-alone  COP  forward  in  the 
windward  plane.  Since  the  load  and  wing  COPs  are 
different  on  the  windward  and  leeward  plane  fins,  this 
results  in  a  net  forward  shift  in  the  COP  for  =  45 
deg  roll  compared  to  the  =  0  deg  computation  of 
Ref.  (1).  This  shift  appears  to  occur  for  all  Mach 
numbers,  is  largest  at  moderate  AOA,  and  goes  to 
zero  at  AOA  of  0  and  90  deg.  At  90  deg,  the 
windward  plane  fins  carry  almost  all  the  load 
compared  to  the  leeward  plane  fins,  but 
geometrically,  the  fins  are  all  aligned  perpendicular 
to  the  AOA  plane. 

Mathematically,  this  geometrical  shift  can  be 
approximated  by:4 


r  + 


(AXcpLb 

cos(3>)  sin(2«) 


V  Cr  +  CtJ  \ 


2 


St 

3 


'  0.8“N 


l  65 


-  0.8 


r  + 


C  +  C. 


;  «  ^  65  (3A) 


\  -  El 

,  2  3 


cos(3>)  sin(2“)  ;  <*  >  65 c 


(3B) 


Here,  CM  was  measured  directly  in  the  data  base  of 
Ref.  9  by  having  the  wing  in  close  proximity  to  the 
body  and  measuring  directly  the  load  on  the  wing  in 
the  presence  of  the  body.  Since  the  normal  force  was 
measured  normal  to  a  single  fin,  to  get  the  normal 
force  on  the  wing  in  the  presence  of  the  body  at 
O  =  45  deg  from  the  data  required  the  data  to  be 
multiplied  by  cos  O.  To  reduce  measurement  errors, 
the  data  from  all  four  fins  were  averaged.  No  attempt 
was  made  to  correct  for  wind  tunnel  errors  near  zero 
AOA  caused  by  flow  misalignments.  These  errors 
can  cause  the  normal  force  curve  to  be  shifted  as 
much  as  a  degree.  This  means  that  the  Cv 

^WfB) 

accuracy  could  have  some  slight  errors  near  zero 
AOA. 


Based  on  the  accuracy  analysis  of  Ref.  1,  fairly 
accurate  values  of  KW(B)  can  be  expected  for  all  but 
the  highest  aspect  ratio  where  the  wing  planform  area 
was  only  about  2  percent  of  the  body  planform  area  in 
the  crossflow  plane.  CN  of  Eq.  (4)  was  arrived  at 
from  Ref  1,  which  in  turn  used  the  data  bases  of 
Refs.  2,  7  and  8. 


The  body-wing  interference  factor  is  defined  as 


K 


B(W) 


NB(W) 


(5) 


Equation  (3)  is  added  to  the  COP  prediction  at 
<I>  =  0  deg1  for  the  roll  orientation  of  45  deg. 

Wing-Bodv  and  Bodv-Wing  Interference  Due  to 
AOA 

The  wing-body  and  body-wing  interference 
factors  were  computed  using  a  combination  of  the 
Ref.  9  data  base  for  0.6  <;  M^  ^  4.6  and  the  Ref.  10 
data  base  at  M  =  0. 1 .  Outside  these  Mach  limits, 
extrapolations  were  made  to  allow  the  methodology 
to  compute  aerodynamics  at  all  Mach  numbers. 

These  extrapolations  were  not  as  difficult  as  they  may 
seem  because  normal  forces  and  COP  have  basically 
leveled  out  at  M  =  4.6,  and  further  increases  in  Mach 
number  produce  fairly  small  changes  in  these 
parameters. 

The  wing-body  interference  factor  is  defined  as 


Unfortunately,  Cw  was  not  a  quantity  that  was 
measured  directly  in  Ref.  7  but  was  computed  from 
three  other  independent  measurements  of  body  alone, 
wing  in  conjunction  with  the  body,  and  total  normal 
force.  The  computation  for  CM  was  then  made  by 

WB(W) 

Cn  =  CN  -  CN  -  CN  (6) 

rtB(W)  W  WB  WW(B)  V  ' 

As  shown  in  the  Ref.  1  error  analysis,  this  process 
gave  potential  errors  that  were  much  higher  than  for 
CN  ,  particularly  for  the  smaller  wings  (AR  ^  1 .0) 
andKigher  Mach  numbers  (M  ^  2.5),  where  the 
CM  term  decreased  to  the  point  where  it  was  within 
the  accuracy  of  the  data.  As  a  result,  much  more 
scatter  in  the  data  is  expected  for  this  term,  and  more 
engineering  judgement  is  required  in  the  empirical 
model  development. 


K, 


W(B) 


Unlike  Ref.  9,  Ref.  10  had  a  fairly  large  wing 
planform  area  compared  to  the  body  planform  area 
(approximately  60  percent).  Moreover,  CN  was 
apparently  measured  separately.  Hence,  CN  could 
be  computed  based  on  direct  measurements,  as  well 
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as  calculated  similar  to  Eq.  (6).  Also,  data  were 
obtained  all  the  way  to  «  =  90  deg  and  at  an  r/s  value 
of  0.25.  As  a  result,  more  confidence  is  placed  on  the 
Ref.  10  body- wing  interference  at  high  AOA  than  the 
Ref.  9  data.  Unfortunately,  the  Ref.  10  data  were 
taken  only  at  MTO  =  0.1,  so  it  is  hard  to  extrapolate  it 
past  about  Mw  =  0.6.  Fortunately,  it  compliments  the 
Mach  number  range  of  the  larger  Ref.  9  data  base 
quite  nicely. 

Figure  2  shows  one  example  of  the  many  figures 
in  Ref.  4  for  the  nonlinear  variations  of  KW(B)  and 
KB(W)  with  AOA  at  M  =  1 .5  for  aspect  ratio  0.5  and  at 
=  45  deg.  Other  Mach  number  and  aspect  ratio 
results  are  given  in  Ref.  4.  It  is  clear  in  examining 
Fig.  2,  and  the  other  figures  in  Ref.  4,  that  inclusion 
of  the  nonlinearities  in  the  interference  factors  is 
essential  in  accurately  developing  a  semiempirical 
nonlinear  aeroprediction  code.  Reference  4  does  this 
similarly  to  Ref.  1,  where  a  mathematical  model  is 
developed  based  on  SBT/LT  plus  deviation  of 
SBT/LT  based  on  data.  That  is, 

^W(B)  ”  [^W(B)^bt  +  (°C» 

^■B(W)  “  [^B(W)]sBT  +  AKfl(W)  (°S  M,  AR,  X) 

LT  V  ' 


The  functions  AKW(B)  and  AKB(W)  are  defined  by  1 1 
tables  in  Ref.  4. 


FIGURE  2.  WING-BODY  AND  BODY- WING 
INTERFERENCE  FACTORS  AS  A 
FUNCTION  OF  AOA 
(M-=  1-5,  r/s  =  0.5) 


Two  of  the  key  issues  associated  with  defining  the 
functions  AKW(B)  and  AKB(W)  are  the  variations  with 
r/s  and  internal  shocks.  Most  of  the  data  available9  is 
for  r/s  of  0.5  with  only  one  set  of  lifting  surfaces 
present.  This  means  any  nonlinearities  associated 
with  r/s  and  internal  shocks  from  forward  to  rearward 
lifting  surfaces  have  to  be  modeled  separately  or  not 
included.  The  most  critical  of  these  issues  appeared 
to  be  the  minimum  value  of  KB(W)  at  high  AOA. 

Figure  3  represents  the  treatment  of  [K^],^  for 
both  the  $  =  0  and  $  =  45  deg  roll  positions.  This 
figure  was  derived  based  on  the  Ref.  9  and  10  data 
bases  for  r/s  =  0.5  and  numerical  experiments  for 
other  r/s  cases.  Referring  to  Fig.  3,  note  that  for  r/s 
^  0.2,  =  0,  there  is  a  difference  in  the  minimum 

value  of  KB(W)  depending  on  whether  there  is  an 
afterbody  present  or  not;  whereas  for  r/s  ^  0.5,  there 
did  not  appear  to  be.  It  should  be  remembered  that 
[Kb(W)]min  is  given  as  a  fraction  of  SBT  in  Fig.  3. 

Thus,  a  value  of  1.0  at  r/s  =  0.2  in  Fig.  3  gives  a  value 
of  of  0.27;  whereas,  a  value  of  0.5  for 

r/s  =  0.5  gives  a  value  of  [K^)]^  of  0.4.  Also 
worthy  of  note  is  the  fact  that  [K^]^  appeared  to 
approach  zero  at  high  Mach  number.  More  data  or 
computational  fluid  dynamics  are  needed  to  more 
precisely  define  Fig.  3. 


0  2  4  6  8 


FIGURE  3.  MINIMUM  VALUE  OF  BODY- WING 
INTERFERENCE  FACTOR  AT  HIGH  AOA 
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Of  equal  importance  in  Fig.  3  is  the  lower  curve 
for  =  45  deg.  Note  that  [KB(W)]iVIIN  is  one-half  to 
one-quarter  of  that  for  <&  =  0  deg.  This  appears  to  be 
the  main  reason  many  missiles  have  lower  normal 
forces  flying  at  O  =  45  deg  than  at  =  0  deg, 
particularly  if  r/s  is  fairly  high  (0.4  or  greater).  Note 
also  the  fact  that  “fin  choking”  is  given  a  separate 
label  for  low  r/s  and  3>  =  45  deg. 

"Fin  choking"  is  a  phenomenon  similar  to  that 
which  occurs  when  an  inlet  becomes  unstarted  or  a 
wind  tunnel  achieves  its  maximum  rate  of  flow  (an 
increase  in  power  produces  no  more  mass  flow 
through  the  inlet).  As  the  body  increases  in  AOA 
with  the  fins  oriented  in  the  “x”  or  cross  orientation, 
the  flow  between  the  fins  will  eventually  "choke"  at 
some  AOA  and  at  moderate  to  large  supersonic  Mach 
numbers.  When  this  happens,  a  strong  shock  is 
formed  just  in  front  of  the  fin,18,19  producing  a  high 
pressure  region  on  the  fins  and  body.  This  high 
pressure  region  is  shifted  forward  from  where  it 
would  be  if  supersonic  flow  occurred  through  the 
fins.  While  the  absolute  value  of  pressure  on  the 
body  is  higher  than  for  the  unchoked  flow,  it  occurs 
over  a  much  smaller  region  and  hence  gives  only 
slightly  higher  body-wing  interference  lift.  Modified 
Newtonian  Theory  was  used  in  Ref.  4  to  define  the 
region  where  “fin  choking”  occurred.  It  was  seen  that 
“fin  choking”  appears  to  occur  on  wing-dominated 
(small  r/s)  configurations  at  3>  =  45  deg  roll,  but  not 
on  body-dominated  configurations  (large  r/s).  A 
fairly  thorough  discussion  of  other  internal  shock 
interactions  was  also  given  in  Ref.  4.  Suffice  it  to  say 
that  additional  work  is  still  needed  before  internal 
shock  interactions  can  be  completely  modeled  with  a 
semiempirical  code. 

Wing-Bodv  and  Bodv-Wing  Interference  Due  to 
Control  Deflection 

The  same  general  approach,  with  slight 
modifications,  to  the  nonlinear  model  of  Ref.  1  for 
0  =  0  deg  roll  is  used  for  the  O  =  45  deg  roll 
position.  In  the  Ref.  1  method,  kW(B)  and  kB(w)  were 
approximated  by 


combined  local  AOA  of  the  wing  |  «*  +  5 1 ,  and  Mach 
number.  The  tables  of  C}  and  C2  for  O  =  45  deg  are 
given  in  Ref.  4. 

In  examining  the  constants  and  model  of  Ref.  4, 
several  physical  phenomena  occur  that  are  modeled  in 
a  semiempirical  sense  by  Eqs.  (8)  and  (9).  These 
phenomena  are  qualitatively  shown  in  Fig.  4.  At  low 
Mach  number,  Fig.  4  indicates  the  SBT  gives  a  low 
value  of  kw(B)  for  small  values  of  «w .  At  a  value  of 
«w  of  about  25  deg,  the  controls  lose  effectiveness  as 
a  result  of  a  combination  of  stall  and  blow-by  effects 
due  to  the  separation  between  the  wing  and  body.  At 
an  «w  of  about  55  deg,  the  controls  have  lost  all 
effectiveness.  At  Mach  numbers  greater  than  about  4, 
the  controls  initially  generate  less  effectiveness  than 
is  generated  by  SBT  for  values  of  «w  up  to  about 
20  or  25  deg.  The  controls  then  become  more 
effective  because  of  nonlinear  compressibility  effects. 
On  the  other  hand,  at  an  «w  of  around  45  to  50  deg, 
the  controls  once  again  begin  to  lose  effectiveness, 
presumably  because  of  shock  interactions  and 
blow-by  effects.  For  Mach  numbers  in  between 
subsonic  and  high  supersonic,  kW(B)  has  behavior  in 
between  the  two  extremes  illustrated  in  Fig.  4. 

CONTROLS  LOSE 
EFFECTIVENESS 
DUE  TO  SHOCK 
INTERACTIONS, 
BLOW  BY  EFFECTS, 
etc 


CONTROLS  LOSE\ 
EFFECTIVENESS  \ 
DUE  TO  STALL,  \ 
BLOW  BY  EFFECTS  \ 

I - 1 - 1 - 1 _ L^K _ I  l| 

0  10  20  30  40  50  60  70 

<*w  (deg) 

FIGURE  4.  QUALITATIVE  TREND  OF  WING- 
BODY  INTERFERENCE  DUE  TO  CONTROL 
DEFLECTION  AS  FUNCTION  OF 
(Mm,  «w) 


CONTROLS 
BECOME  MORE 
EFFECTIVE  DUE  TO 
COMPRESSIBILITY 


kw(B)  C\(M)  [kW(B)^BT  +  C2(  |  ocw  | ,  M)  (8) 
kB(W)  =  [kB(W)^BT  (9) 

The  parameters  Cj  and  C2  were  derived  based  on 
numerical  experiments  of  the  AP95  compared  to  data. 
This  was  because  many  of  the  fins  in  the  Ref.  9  data 
base  were  too  small  to  allow  accurate  estimates  of 
kw(B)  as  a  function  of  parameters  of  interest.  As  a 
result,  total  missile  load  data  were  used  and  empirical 
values  of  C}  and  C2  estimated  as  a  function  of 


In  comparing  the  nonlinear  control  deflection 
models  for  $  =  0  and  45  deg  roll  in  Ref.  4,  a  lot  of 
similarity  is  seen.  The  constants  for  $  =  45  deg  are 
slightly  different  than  those  for  O  =  0  deg  and  the 
values  of  «w  where  the  nonlinearities  begin  are 
somewhat  different.  However,  by  and  large,  Eq.  (8) 
holds  for  both  the  $  =  0  and  45  deg  roll  cases.  It 
should  be  pointed  out  that  in  Ref.  1,  mostly  linear 
variations  of  kw(B)  with  «w  were  used.  However,  these 
have  been  improved  upon  for  the  3>  =  45  deg  case 
with  cubic  fits  of  control  deflection  data.  As  such,  all 
nonlinear  effects  are  included  in  the  variations  of 
kW(B)  as  a  function  of  Mach  number  and  | «  +  8 1 .  It 
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should  also  be  noted  that  kw(B)  and  kB(w)  are  multiplied 
by  1.414  to  indicate  that  all  four  fins  are  assumed  to 
be  deflected  by  an  equal  amount  in  the  =  45  deg 
roll  position. 

Reference  4  also  derived  new  semiempirical 
models  for  wing-tail  interference.  This  is  the  subject 
of  another  paper  and  is  thus  not  discussed  here. 
However,  the  comparisons  with  data  in  the  Results 
and  Discussion  section  have  the  new  wing-tail 
interference  model  included  when  the  configuration 
has  two  sets  of  lifting  surfaces  present. 

Results  and  Discussion 

The  new  nonlinear  aeroprediction  methodology 
for  the  4>  =  45  deg  roll  position  has  been  validated 
against  many  configurations  within  and  outside  the 
data  bases  upon  which  the  methodology  was 
developed.  Reference  4  gives  results  for  seven  wing- 
or  tail-body  cases  and  five  wing-  or  canard-body-tail 
cases.  Three  cases  are  chosen  here  to  illustrate  the 
accuracy  of  the  new  methodology  as  well  as  illustrate 
remaining  problems  in  the  state-of-the-art  in  both 
experimental  testing  and  semiempirical  code 
development.  Two  of  these  cases  are  given  in  Ref.  4 
in  more  detail,  and  the  third  case  has  not  been  shown 
previously  for  the  <I>  =  45  deg  roll  aerodynamics. 

The  first  of  these  cases  is  from  Ref.  16  and  is 
referred  to  as  the  air-slew-demonstrator-vehicle.  Data 
were  available  at  both  <£  =  0  and  45  deg  roll  at 
M  =  0.6  to  1.3  and  to  AOA  90  deg.  Figure  5  shows 
the  configuration  schematically  and  the  normal  force 
coefficient  results  for  M  =  0.6,  0.9,  and  1.3.  This 
case  was  chosen  for  several  reasons.  First,  it 
illustrates  the  accuracy  of  the  code  at  both  $>  =  0  and 
45  deg  at  moderate  subsonic  Mach  numbers  to  AOA 
90  deg.  Secondly,  it  illustrates  a  problem 
encountered  in  validating  the  code  at  subsonic  and 
transonic  Mach  numbers  based  on  wind  tunnel  data. 
This  problem  is  wind  tunnel  model  sting  interference. 

Reviewing  some  of  the  wind  tunnel  model  support 
interference  literature  (Refs.  17-20),  several 
conclusions  were  reached.  These  were  that:  a)  For 
low  Mach  number,  the  problem  of  estimating 
interference  effects  of  a  strut  or  sting  mount  on  the 
model  aerodynamics  at  high  AOA  are  not  known 
precisely — this  is  still  true  today,  b)  The  preferable 
mount  between  a  sting  and  a  strut  at  high  AOA  is  the 
sting,  c)  The  sting  tends  to  give  positive  interference 
(CN  too  high)  and  the  strut  negative  interference  (CN 
too  low). 

The  affected  region  seems  to  be  in  the  AOA  range 
30-80  deg.  Sting  CN  values  can  be  high  as  much  as 
10-15  percent,  and  strut  CN  values  for  struts  mounted 
in  the  mid-body  region  can  be  low  by  as  much  as 


25-30  percent.  For  these  reasons,  the  Improved 
Aeroprediction  Code  1995  (IAP95)  methodology  is 
intentionally  designed  to  underpredict  normal  force 
on  configurations  at  subsonic  Mach  numbers  at  high 
AOA,  where  test  data  are  from  a  sting  mount  and 
therefore  it  is  also  expected  the  predictions  will  be 
higher  than  those  from  a  strut-mounted  model.  The 
under/overprediction  problem  on  CN  at  high  « 
appears  to  go  away  at  Mach  numbers  slightly  greater 
than  one.  This  is  suspected  to  be  due  to  the  reduction 
in  the  upstream  influence  of  the  sting  on  the  body  and 
the  reduction  of  the  wake  effect  on  the  body  (strut 
mount). 

Since  the  Fig.  5  configuration  was  mounted  in  the 
tunnel  with  a  strut,  this  is  suspected  to  be  the  major 
reason  for  the  overprediction  of  the  IAP95 
methodology  compared  to  data  at  M  =  0.6,  0.9  and 
AOA  >  50  deg.  While  COP  predictions  are  not 
shown,  results  are  within  the  ±  4  percent  of  body 
length  goal  at  all  AOA  and  Mach  numbers.  This 
implies  any  loss  of  normal  force  due  to  the  strut  is 
distributed  all  along  the  body.  Also,  if  one  excludes 
the  AOA  >  50  deg  comparison  at  M  =  0.6  and  0.9 
where  wind  tunnel  measurement  errors  are  in 
question,  the  IAP95  is  well  within  the  ±10  percent 
accuracy  goal  on  normal  force  as  well. 

The  last  two  configurations  considered  in  the 
validation  process  utilize  all  the  new  nonlinear 
methodology  developed  for  both  the  3>  =  0  and 
45  deg  roll  orientations.  The  first  of  these  is  a 
version  of  the  SEASPARROW  missile  as  shown  in 
Fig.  6,  with  wind  tunnel  results  taken  from  Ref.  21 . 
This  configuration  has  fairly  large  wings  and  tails, 
with  the  wings  used  as  the  control.  Data  were  taken 
at  Mach  numbers  of  1.5, 2.0,  2.35, 2.87,  3.95  and 
4.63.  Figs.  6B,  6C,  and  6D  give  comparisons  of  the 
IAP95  with  the  data  at  M  =  1.5,  2.87,  and  4.63 
respectively.  Results  are  shown  for  control 
deflections  of  10  deg  for  M  =  1 .5  and  20  deg  for 
M  =  2.87  and  4.63.  Also,  all  results  for  both  the 
€>  =  0  and  45  deg  roll  positions  are  shown.  Several 
comments  are  in  order  with  respect  to  the  overall 
comparisons.  First  of  all,  the  LAP95  model  achieves 
its  goal  of  predicting  average  accuracy  of  CA,  CN  of 
±10  percent  and  XCP  of  ±4  percent  (?B  on  this 
configuration  at  both  the  3>  =  0  and  45  deg  roll 
positions.  Secondly,  CN  is  predicted  equally  well  at 
4>  =  0  and  45  deg.  However,  CM  is  predicted  better  at 
$  =  45  deg  than  at  =  0  deg  due  to  the  COP  shift 
discussed  earlier.  No  such  shift  has  been  applied  at 
$>  =  0  deg.  Apparently  one  is  needed,  but  it  is  not 
clear  what  the  physical  justification  is.  As  seen  in  the 
pitching  moment  predictions  for  the  4>  =  0  deg  roll 
orientation,  the  IAP95  in  general  is  slightly  too  stable. 
The  third  point  is  that  CA  prediction  is  slightly  better 
for  the  <E>  =  0  deg  roll  orientation  than  the  3>  =  45  deg 
position.  Apparently,  the  factor  of  1.414  applied  to 
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FIGURE  5.  NORMAL  FORCE  COMPARISONS  OF  THEORY  AND  EXPERIMENT 
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the  fins  for  the  control  deflection  component  of  axial 
force  is  quite  adequate  at  low  AOA  at  all  Mach 
numbers,  as  seen  in  the  figures.  However,  at  high 
AOA,  this  factor  appears  to  be  too  high  at  the  lower 
supersonic  Mach  numbers  and  too  low  at  the  higher 
supersonic  Mach  numbers.  It  is  suspected  that  at  the 
lower  supersonic  Mach  numbers,  only  the  fins  in  the 
windward  plane  should  have  the  full  factor  of  yfl, 
whereas  those  in  the  leeward  plane  should  have  a 
lower  factor.  For  high  supersonic  conditions,  it  is 
suspected  the  bow  shock  and  internal  shock 
interactions  actually  add  to  the  factor  of  yjl .  An 
empirical  model  for  the  control  deflection  component 
of  CA  to  account  for  this  physics  would  improve  the 
CA  comparison,  but  time  did  not  permit  this  effort. 


FIGURE  6A.  AIR-TO-AIR  MISSILE 

CONFIGURATION  USED  IN  VALIDATION 
(FROM  REF.  21) 

The  final  point  to  be  made  is  concerning  the  high 
Mach  number,  high  AOA  conditions  in  Fig.  6. 

Notice  that  the  pitching  moment  suddenly  loses 
stability  above  AOA  of  30  deg  and  <$  =  0  deg  roll  at 
M  =  4.63.  The  normal  force  also  decreases  somewhat 
as  well.  It  is  believed  this  is  due  to  internal  shock 
waves  from  the  bow  shock  and  wing  shock 
intersecting  the  tail,  causing  a  loss  of  tail  lift  and  a 
sudden  decrease  in  stability.  No  accounting  of  these 
internal  shocks  from  a  forward  mounted  fin  to  an  aft 
mounted  fin  is  made  in  the  IAP95.  It  is  also 
interesting  to  note  that  above  AOA  70  deg,  while  the 
results  are  not  shown,  the  IAP95  predictions  for 
pitching  moment  agree  quite  well  with  the  data. 
Apparently  the  internal  shock  interactions  are 
important  between  AOA  of  about  30  to  60  deg. 

The  last  configuration  chosen  for  validation  of  the 
new  methodology  against  wind  tunnel  results  is 
shown  in  Fig.  7A.  The  data  as  well  as  DATCOM3 
results  are  taken  from  Ref.  22.  The  control  is  from 
the  canards.  The  configuration  is  over  22  cal  in 
length  with  aspect  ratio  tails  of  0.9  and  canards  1.57. 
Data  were  taken  for  Mach  number  of  0.2  but  to  AOA 
of  50  deg  with  control  deflections  of  0  and  ±20  deg. 
The  ±20  deg  control  deflection  cases  are  illustrated  in 
Fig.  7B.  Both  the  IAP95  and  Missile  DATCOM 
results  from  Ref.  22  are  shown  along  with  the  data  for 


M=1.5.  5  =  10° 


FIGURE  6B.  AXIAL,  NORMAL,  AND 
PITCHING  MOMENT  COEFFICIENT 
COMPARISON  OF  THEORY  AND 
EXPERIMENT 
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M=4.63.  8  =  20° 


M=2.87.  5  =  20° 


FIGURE  6C.  AXIAL,  NORMAL,  AND 
PITCHING  MOMENT  COEFFICIENT 
COMPARISON  OF  THEORY  AND 
EXPERIMENT 


FIGURE  6D.  AXIAL,  NORMAL,  AND 
PITCHING  MOMENT  COEFFICIENT 
COMPARISON  OF  THEORY  AND 
EXPERIMENT 
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<5  =  45  deg  roll  only.  In  general,  the  IAP95  gives 
very  good  comparisons  with  data  and  appears  to 
exhibit  most  of  the  nonlinearities  of  the  data.  The 
6  =  +20  deg  IAP95  predictions  agree  with  the  data 
slightly  better  than  5  =  -20  deg.  Missile  DATCOM 
results  also  give  reasonably  good  agreement  with 
data.  Reference  1  compared  the  AP95,  Missile 
DATCOM  and  data  for  this  same  case  at  <5  =  0  deg 
roll.  The  AP95  results  were  about  as  good  at  <5  =  0 
as  <5  =  45.  However,  the  Missile  DATCOM  gave 
better  results  compared  to  data  at  <5  =  45  deg  versus 
<5  =  0  deg.  It  is  not  known  whether  this  statement 
holds  true  in  general,  however. 


T 

5.0 

FIGURE  7A.  CANARD  CONTROLLED  MISSILE 

CONFIGURATION  (ALL  DIMENSIONS  IN 
INCHES,  FULL  SCALE) 

Summary.  Conclusions  and  Recommendations 

New  technology  has  been  developed  to  allow 
engineering  estimates  of  aerodynamics  of  most 
tactical  weapon  concepts  for  the  <5  =  45  deg  roll 
position  (fins  oriented  in  "x"  or  cross-fin 
arrangement).  New  technology  developed  for  the 
<5  =  45  deg  roll  position  and  discussed  in  this  paper 
includes: 

a)  Nonlinear  wing-body  and  body-wing 
interference  factor  methodology  due  to  AOA 
and  control  deflection. 

b)  An  approximate  method  to  estimate  wing- 
alone  shift  in  COP  at  <5  =  45  deg  and  high 
AOA. 

c)  An  approximate  method  to  account  for  "fin 
choking"  at  high  Mach  number  and  AOA 

This  new  technology,  along  with  the  <5  =  0 
methodology  of  Ref.  1,  now  allows  trim 
aerodynamics  to  be  computed  in  the  two  roll- 
stabilized  planes  that  most  of  the  world’s  missiles  fly. 
The  data  bases  upon  which  the  methodology  was 
based  were  limited  in  Mach  number  to  4.5  and  AOA 
of  40  deg.  However,  engineering  judgement  and 
other  data  were  used  to  allow  calculations  to  be 
performed  to  AOA  90  deg,  Mach  numbers  up  to  20, 
for  axisymmetric  solid  rocket  weapons  with  up  to  two 
sets  of  lifting  surfaces. 


Based  on  the  new  methodology  and  computations 
to  date,  the  following  conclusions  are  made. 

a)  The  two  primary  reasons  missile  normal 
forces  are  lower  at  <5  =  45  deg  than  at 

<5  =  0  deg  are  that  the  minimum  value  of 
body  carryover  lift  at  high  AOA  is  lower  for 
<5  =  45  than  <5  =  0.  The  larger  the  value  of 
r/s,  the  more  difference  between  these 
minimum  values.  If  the  configuration  has 
two  sets  of  lifting  surfaces,  the  wing-tail 
interference  is  higher  for  AOA  >  20  deg  for 
the  <5  =  45  deg  plane  than  for  the  <5  =  0  deg 
plane,  also  contributing  to  a  lower  normal 
force  for  <5  =  45  deg  and  a  less  stable 
configuration. 

b)  At  high  AOA:  KW(B)  approaches  1.0  for 
most  Mach  numbers;  KB(W)  approaches  some 
minimum  value  that  is  a  function  of  r/s  and 
Mach  number;  kw(B)  and  kB(w)  are  nonlinear 
in  total  local  AOA  on  the  wing  and  are 
functions  of  Mach  number  and  total  AOA  on 
the  wing. 

c)  "Fin  choking"  appears  to  be  more  of  a 
problem  on  larger  fin  (small  r/s) 
configurations  at  O  =  45  deg  than  on  smaller 
fin  cases  (large  r/s). 

d)  For  Mach  numbers  less  than  one  and  AOA  > 
30  deg,  it  is  not  clear  what  the  correct  values 
of  experimental  normal  force  are  for  a  given 
configuration,  based  on  wind  tunnel  results 
available  in  the  literature. 

e)  Internal  shock  interactions  between  forward- 
mounted  and  aft-mounted  fins  become 
increasingly  important  as  both  AOA  and 
Mach  number  increase.  The  current 
methodology  does  not  account  for  these 
effects. 

f)  In  general,  ±10  percent  average  accuracy  has 
been  maintained  for  both  normal  and  axial 
force  coefficients  and  ±4  percent  of  body 
length  for  COP  in  the  <5  =  45  deg  roll 
position.  Exceptions  to  this  are  at  subsonic 
Mach  number  and  high  AOA  where  data 
accuracy  is  in  question,  and  at  high  Mach 
number  and  high  AOA  for  configurations 
that  have  two  sets  of  lifting  surfaces,  where 
internal  shock  interactions  may  be  a 
problem. 

g)  The  current  overall  approach  of  using  linear 
theory,  slender  body  theory,  or  second  order 
theory  for  low  AOA  aerodynamics,  and 
estimating  the  nonlinear  aerodynamic  terms 
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FIGURE  IB.  AXIAL,  NORMAL,  AND  PITCHING  MOMENT  COEFFICIENT 
COMPARISONS  OF  THEORY  AND  EXPERIMENT 
(M  =  0.2,  *  =  45  DEG) 
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individually  and  directly  from  wind  tunnel 
data  bases,  appears  to  be  the  key  to  the 
above-mentioned  average  accuracy  levels. 

Based  on  this  and  previous  research  (Refs.  1,  6), 
the  following  recommendations  are  made  for 
additional  work: 

a)  Additional  wind  tunnel  measurements  or 
computational  fluid  dynamics  cases  need  to 
be  made  to  help  define  nonlinearities  of  the 
interference  effects  as  a  function  of  r/s. 

b)  A  method  is  needed  to  accurately  correct 
wind  tunnel  data  at  subsonic  Mach  number 
and  high  AOA  for  sting  or  strut  mounting 
effects. 

c)  An  engineering  method  is  needed  to  estimate 
internal  shock  interaction  effects  for  high 
AOA  and  Mach  number. 

d)  Any  future  wind  tunnel  test  for  measuring 
component  aerodynamics  should  be  done 
with  lifting  surfaces  large  enough  to  separate 
out  body  and  wing  lift  accurately,  with  wings 
mounted  in  the  middle  of  the  body  and, 
preferably,  with  simultaneous  measurements 
of  body  forces  in  conjunction  with  the  wing 
and  wing  forces  in  conjunction  with  the 
body.  This  would  allow  more  accurate 
determination  of  the  interference  terms 
directly,  without  subtraction  of  two  large 
numbers  to  obtain  a  small  term. 
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